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Abstract 

Bound state properties of few single and double-A hypernuclei is critically ex¬ 
amined in the framework of core-A and core-l-A -|- A few-body model applying 
hyperspherical harmonics expansion method (HHEM). The AA potential is chosen 
phenomenologically while the core-A potential is obtained by folding a phenomeno¬ 
logical AN interaction into the density distribution of the core. The depth of the 
effective AN potential is adjusted to reproduce the experimental data for the core- 
A subsystem. The three-body Schrodinger equation is solved by hyperspherical 
adiabatic approximation (HAA) to get the ground state energy and wave function. 

The ground state wavefunction is used to construct the supersymmetric partner 
potential following prescription of supersymmetric quantum mechanics (SSQM) 
algebra. The newly constructed supersymmetric partner potential is used to solve 
the three-body Schrodinger equation to get the energy and wavefunction for the 
first excited state of the original potential. The method is repeated to predict en¬ 
ergy and wavefunction of the next higher excited states. The possible number of 
bound states is found to increase with the increase in mass of the core of the hy¬ 
pernuclei. The Root Mean Squared (RMS) matter radius and some other relevant 
geometrical observables are also predicted. 
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I Introduction 


Since the discovery of the light exotic hypernuclei in the early sixties [1-2], peoples are 
giving much attention to study the structure of such multistrange hypernuclei [3-4]. Some 
of the experimentally observed one- and two-A hypernuclei are ^He, AAHe, iBe, AA^e, 
JfC and [5-11]- Information about the AA and AN forces can be extracted from the 
bound state properties of such nuclei. These informations on AA and AN interaction to¬ 
gether with the well known NN interaction may give better understanding of the strong 
quark-quark (qq) interaction, since the hyperons also have qqq structure as that of nu¬ 
cleons. As for example proton has three quarks namely up,up,down (p —)■ uud), neutron 
has quarks up,down,down (n —)■ udd), A° has quarks up,down,strange (A° —)■ uds) etc. 
In the earlier stage of their discovery, emulsion experiments provided some information 
on the binding energies of A-particle in the light exotic hypernuclei and their decay rates 
(life times~ 10“^°sec) [2]. From the available binding energy data physicists gathered 
some qualitative informations about the A-nucleon (AA^) interaction and single particle 
potential strength for the A-particle in hypernuclei [12]. The hyperon-nucleon scattering 
experiments have also been performed but these are not so profound and are still in the 
primary stages and do not give detailed phase shifts to construct the potential reliably. 
Some An and T,N total cross-sections and very few angular distribution at low energies 
have been measured [13-18], but these are not sufficient to allow the phase shift analysis. 
Under the circumstances, the bound state properties of single- and double- A hypernuclei 
can only give useful indirect information about AN and AA interactions. For example, 
one can take phenomenological forms of AN and AA interactions and see if they repro¬ 
duce the experimentally observed observables of the hypernuclei. Alternatively one can 
adjust the parameters of the empirical potential to reproduce the bound state properties 
and thus predict the effective AN and AA interactions. In some earlier works [19-22] 
physicists have used variational and approximate few body methods for the hypernuclei 
treating them a few body system. In the present work, we investigate the structure and 
properties of the ground and excited states of double-A hypernuclei with mass number 
A = 6,8,34,42,92 and ground state of their binary subsystems consisting the core and 
one outer core A hyperon. A simple two-body model has been adopted here for the single- 
A hypernucleus {A = 1) consisting of a Ac = 4, 6, 30,40, 90 core and one valence A 

particle. And for the double-A hypernuclei (A = Ac + 2), we adopt a three-body model 
consisting core of mass Ac (=4,6,30,40,90) and two valence A particles. For AA potential 
a phenomenological form has been chosen. The core-A potential has been obtained by 
folding an effective AN potential by the density of the core nucleus. The strength of 
the effective AN potential for a given core is adjusted to reproduce the experimentally 
known BE of the single A hypernuclei |He, ^He [6], ]^Be [23], [24], ^^Ca [24] and ®^Zr 

[25]. The same folded core-A potential is then used for the double-A hypernucleus for 
the same core. 

We used hyperspherical harmonics expansion (HHE) method to solve the three-body sys¬ 
tem. This method is a powerful tool for the ab initio solution of the few body Schrodinger 
equation for a given set of binary interaction potentials among the constituent particles. 
This method has been used for bound states in atomic [26-33], nuclear [34-40] and par¬ 
ticle physics [41-43]. Attempts have been made to use it in scattering problems as well 
[44]. In this method, the wave function is expanded in a complete set of hyperspheri¬ 
cal harmonics (HH), which are, for a three-body system, the six-dimensional analogue 
of ordinary spherical harmonics. The Schrodinger equation reduces to a set of coupled 
differential equations which can be solved numerically. The HHE method is essentially 
exact and involves no approximation other than an eventual truncation of the expan¬ 
sion basis. Any desired precision in the binding energy can, in principle, be achieved by 
gradually expanding the expansion basis and checking the rate of convergence. However 
the number of coupled differential equations and hence, the complexity in the numerical 
solution increases rapidly as the expansion basis is increased by introducing larger hyper 
angular momentum quantum numbers. Computer limitations set an ultimate limit to 
the precision attainable. 
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We have calculated the two-A separation energy (-Baa), AA bond energy (A.Baa) and 
some size parameters for all the above double A hypernuclei. Here .Baa and AA bond 
energy A.Baa are dehned as 

BU^^Z) = [M(^-^Z) + 2AA - M(i^Z)]c-^ (1) 

and 

AHaa = 5aa(^a^) - 25a(^'Z) (2) 

As the ground state BE is quite large, one may expect the excited states to be observable 
experimentally, although they have not been reported so far. However, the HHE method 
is known for its rather slow convergence, which is especially manifested in the excited 
states. The excited state wave function is more spread out than the ground state and 
therefore needs larger Kmax for the same degree of convergence. We circumvent this 
difficulty by invoking supersymmetric quantum mechanics (SSQM) [45]. In SSQM, we 
construct the supersymmetric partner potential (V 2 ) of the original potential (W) [1]. 
According to SSQM, the spectra of Vi and V 2 are identical except for the absence of the 
ground state of Vi in the spectrnm of V 2 [45]. Hence the gronnd state of V 2 corresponds 
to the hrst excited state of Vi. Thns we can solve for the gronnd state of V 2 , whose 
convergence is fast (being a gronnd state). Its BE is the BE of the first excited state of 
the original potential. The wave fnnction of the hrst excited state of Vi can be obtained 
by applying an appropriate operator on the gronnd state wave fnnction of V 2 . 

The paper is organized as follows: In section H, we review the HHE method for a three- 
body system consisting of non identical particles with a brief review of SSQM. Results 
of calculation and discussion are presented in section HI. Finally in section IV we draw 
our conclusions. 


II Use of hyperspher spherical harmonics expansion 
method and super symmetric quantum mechanics 


We treat each of the six double-A hypernnclei ivHe, i^Be, AA^e, aaS, AA^a, AA^r as a 
three body system, where each of the core ^He, ®Be, ®He, ^^S, ^'^Ca, ®°Zr is labeled as 
particle 1 and the two valence A hyperons as particles 2 and 3 respectively (see Fig. 1). 

For a given partition i (in which i is the spectator), a set of Jacobi co-ordinates is 
dehned as: 


Xi = 

rtijirikM 

1 

4 


- n) 


2 

Vi = 

rrijinj^M 

4 

(a 

mjfj+mkr‘k\ > 
rrij+rrik ) 


-R = A + mjrj + rrikn). 


( 3 ) 


for (i,j,k) = (l,2,3 cyclic). Here rrii, fl are the mass and position of the particle and 

M = mi + mj + mk- R is the coordinate of the centre of mass of the system. Since the 
interactions among the constituent particles depend only on their relative separations, 
the centre of mass motion separates ont antomatically. Thns the relative motion of the 
three-body system is described by the Schrodinger eqnation 


-^(V^, + Vjj + Vjk{xi) + Vki{xi, Vi) + Vij{xi, Vi) - E 


'^{xi.Vi) = 0 (4) 
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is an ehective mass parameter and W is the interaction potential 
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between particles i and j. In terms of the hyperspherical variables [40] 
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Figure 1: Choice of Jacobi coordinates for the partition 1. 


the Schrodinger equation becomes 


1 d ^ d 
2 p ^ ^ 


/c2(a) 


} + V{p,n,)-E 


VI/(p, a) 


0 , 


( 6 ) 


where the hyper-radius p is invariant under three-dimensional rotations and permutations 
of the particle indices. The hve hyper angles constituted by polar angles of 

Xi), polar angles {0yi,4>yi) of pi and the angle <Fj, are collectively denoted by —)■ 
6 * 3 ;., 0J.., 0yj}. Thus (p, fJj) constitutes the six hyperspherical variables of which 

the hve hyper angles Qi depend on the choice of partition i. In the above equation 
V (p, ffj) = Vjk{xi) + Vki{xi, yl) -|- Vij{xi, Pi) is the total interaction potential expressed in 

terms of the hyperspherical variables in the partition i and is the square of hyper 

angular momentum operator given by [40] 






4 cot 


d 




(7) 


where and are the squares of ordinary orbital angular momentum operators 

associated with Xi and yl motions respectively. The operator /C^ satishes an eigenvalue 
equation [40] 

= K{K + 4)yKaA^^). (8) 

Here K (= 2ni + + lyi,ni being a non-negative integer) represents the hyper angular 

momentum quantum number and a* = {lx^, ly^, L, M} is a short hand notation, in which 
L and M denote the total orbital angular momentum and its projection, respectively. 
The normalized eigen functions yxaii^i) sxe called the hyperspherical harmonics (HH) 
and these depend on the choice of partition. A detailed analytic expression for yRaA^i) 
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can be found in [40]. In the HHE method, \I/(p, hlj) is expanded in the complete set of 
HH corresponding to the partition i: 




Kai 


P 


( 9 ) 


Substitution of Eq. (9) in Eq. (6) and the use of orthonormality of HH leads to a set of 
coupled differential equations (CDE) in p 


_h^ (jp_ _ 
2^ \ dp^ 




-E 


UKaiip) 


+ EiC'a,' < Kai I V{p,Vti) I K'a- > UK'a/{p) = 0, 

where Ck = K + and 


( 10 ) 


< A'a,|l'(p,Sl,)|A"Q/>= / (11) 

JUi 

The number of coupled differential equations to be solved, Eq. (10) is in principle, an 
inhnite one, which arises out of an inhnite number of basis states Eq. (9). For practical 
computation, the hyperspherical harmonics expansion basis in Eq. (9) is truncated to 
a hnite set by retaining all values of K up to a maximum Kmax- For a given K, all 
allowed values of the set of quantum numbers a* are included. The basis states is further 
reduced by imposing the constraints arising out of symmetry requirements and associated 
conserved quantum numbers. Thus the infinite set of CDE is now reduced to a finite set. 
This hnite set of CDE is solved numerically by hyperspherical adiabatic approximation 
[31,46] which is based on the idea of adiabatic separation of hyper-angular motions from 
the hyper-radial one. In this procedure the potential matrix < Kai\V\K'a[ > together 
with the diagonal hyper-centrifugal term is diagonalized for a hxed value of p. The lowest 
eigenvalue Ci;o(p) and the corresponding eigenvector XKaipip) are obtained as parametric 
function of p. Then the set of CDE, Eq. (10) is approximately decoupled into a single 
ordinary differential equation 


'_h^ (P 

2p dp"^ 


+ ujq{p) - E 


P{.p) = 0 


( 12 ) 


Numerical solution of Eq. (12), subject to appropriate boundary conditions gives the 
ground state energy E and partial waves of the three-body system. The partial waves 
are given by 

Ukc^Xp) - Po{p)XKaiM ( 13 ) 

We now present a brief review of properties of supersymmetric quantum mechanics 
(SSQM) [45] used here to calculate the observables for the excited states. Let us consider 
a one dimensional Schrodinger equation of the form 


2/i dp"^ 


+ V{p) - E 


P{.p) = 0 


(14) 


For the single-A hypernuclei, V(p) is just the A-core potential plus the centrifugal term 
and p is the core-A separation. For double-A hypernuclei V(p) is the lowest eigen potential 
LJo(p) and p is the hyper-radius. Suppose the ground state energy is Eq. We shift the 
energy scale by Eq, such that the ground state is at zero energy and the corresponding 
potential is renamed Vpp) {=V{p) — Eq). Then for the ground state 


2p dp‘^ 


+ r,(p) 


i’oip) = 0 


(15) 
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From this we get 



' ^ 2pMp} 

( 16 ) 

Let us now dehne a 

super potential, W (p) as 



Wip)^ 

(17) 

and two operators 

^ + \ 

= -^,i + W{p) / 

( 18 ) 

Then one can easily check that 



VM = w\p) - -^w’(p) 

( 19 ) 

and 

+ 

T, 
^ hs 

1 

III 

(20) 

Here Hi is the original Hamiltonian in the shifted energy scale. 
Hamiltonian 

We next dehne a partner 

(21) 


in which the partner potential is given by 


V,(p) = W\p) + ( 22 ) 

Let us denote the n— eigenvalue and corresponding eigen function of H 2 by and 
ipn\p) respectively, while the n— eigenvalue and corresponding eigenfunction of Hi are 
denoted by En and t/’n(p) respectively. Then it is easily verihed that 

= AA^AiJr. = AHi^j^ = E^A^;^), (23) 

which means that Aipn is an eigenfunction of H 2 corresponding to energy En- Similarly 

= E^{A^i>^;^) ( 24 ) 

Thus is an eigenfunction of Hi corresponding to energy E^^K Now from Eq. (17) 

and (18) 

71^70 = 0. (25) 

Thus AipQ is a trivial solution of Eq. (23). The ground state of H 2 will thus be Tl'^i 
corresponding to energy Ei (note that Eq < Ei < E 2 < Eq < .). Hence Eq^^ = Ei, 

and from Eq. (24), '^icx)Hi'^o^^.Thus if we solve for the ground state of V 2 (p), we get the 
energy Ei (in shifted energy scale), which is the energy of the hrst excited state of Vi{p). 
It is easily seen that the corresponding normalized eigenfunction is 

^i(p) = 

Hence we numerically solve Eq. (14) for the ground state and calculate the super¬ 
potential W{p) and the partner potential V 2 (p) according to Eq. (17) and (22) respec¬ 
tively. Then we solve the same Schrodinger equation with V 2 (p) to obtain E^'^ and 
From ip^\ we get t/’i(p) using eq(26). Finally the energy of the hrst excited state of Vi{p) 
is obtained by back shifting Eq . 
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Ill Results and Discussions 


In the present calculation we have taken the core to be structureless. Since the core 
(^He, ®Be, ®He, ^“^Ca, ^‘^Zr) contains only nucleons and no A-particles, there is no 
symmetry requirements under exchange of the valence A particles with the core nucleons. 
The only symmetry requirements are (i) anti-symmetrization of the core wave function 
under exchange of the nucleons and (ii) anti-symmetrization of the three body wave 
function under exchange of the two A particles. The former is implicitly taken care of 
by the choice of the core as a building block. For double-A hypernuclei, the latter is 
correctly incorporated by restricting the values, explained in the following. Thus, 
within the three-body model, the symmetry requirements are correctly satished without 
any approximation. For the single-A hypernuclei, there are no symmetry requirements. 
The ground state of all experimentally known double-A hypernuclei have a total angular 
momentum quantum number J = 0 and even (or positive) parity. We assume this to be 
true for all double-A hypernuclei with core having N = even, Z = even. The possible 
total spin (S) of the three-body system (core -|- A -|- A) can take two values 0 or 1 since 
the spin of the core in all the above cases has a value 0. Thus the total orbital angular 
momentum L can be either 0 or 1 corresponding to S=0 or 1 respectively. Hence the set of 
quantum numbers (LS)J for the ground state of all even-even core double-A hypernuclei 
is (00)0 and (11)0 which corresponds to and states. Since the ground state of 
these nuclei have dehnite parity, odd and even parity states will not mix. In our case, 
the ground state of all the six double-A hypernuclei is a pure state. Since the core 
is spin less, the spin singlet state (S=0) corresponds to zero total spin of the valence 
A- particles (i.e. >S'23=0). Hence the spin part of the wave function is antisymmetric 
under the exchange of the spins of the two A-particles. Thus the spatial part must be 
symmetric under the exchange of the two A-hyperons. The symmetry of the spatial part 
is determined by the hyper spherical harmonics, since the hyper radius p and hence the 
hyper radial partial waves {Uxaip)) are invariant under permutation of the particles. 
Under the pair exchange operator P 23 which interchanges particles 2 and 3, aq —)■ —Xi 
and yl remains unchanged (see Eq. (3)). Consequently P 23 acts like the parity operator 
for (23) pair only. Choosing the two valence A-hyperons to be in spin singlet state 
(spin antisymmetric), the space wave function must be symmetric under P 23 . This then 
requires Ixi to be even. For the spin singlet state total orbital angular momentum, L=0, 
hence we must have = ly^ = even integer. Since K = 2ni + where ni is a 

non-negative integer, K must be even and 

4i = /j/i = 0,2,4,.., P/2 ifK/2iseven \ 

0,2,4,.., (P/2 - 1) if K/2 is odd j ■ 


Again for the triplet state (S=l), the two valence A-hyperons will be in spin triplet state 
(*S '23 = 1, spin symmetric). Hence the space wave function must be antisymmetric under 
P23. This then requires to be odd. For the spin triplet state, the total orbital angular 
momentum, L=l, hence may take values Ix^ and Uj ± 1 (and only 1 if Ix^ = 0) but the 
parity conservation allows lyi=lxi only (except Ix^ = 0). Again since P = 2ni + Ix^ -|- ly^, 
where ni is a non-negative integer, K must be even (P 7 ^= 0) and 


Ixi = lyi = 1,3,5,.., P/2 if K/2 is odd 

1, 3,5,.., (P/2 — 1) if P/2 is even 


(28) 


For a practical calculation, the HH expansion basis (eq(9)) is truncated to a maximum 
value {Kmax) of K. For each allowed K < K^ax with K=even integers, all allowed values 
of Ixi are included. The even values of correspond to L=0, S=0 and odd values of 
Ixi correspond to L=l, S=l. This truncates eq(lO) to a set of N coupled differential 
equations, where 

N = f + 1 ) (A^ + 1 ) if Kraaxl‘^ IS oven 1 

+ 2 ) if K^ax/‘2 is odd j ’ 
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which will be solved by the hyperspherical adiabatic approximation (HAA) [46]. 


Since realistic AA potential is not available at this stage. We used phenomenological three 
term Gaussian AA potential with parameters adjusted to reproduce the experimental data 
for known double-A hypernuclei. The AA potential used here is given by 


3 

bAA(r) = Xl^iexp(- 

i=\ 



(30) 


The parameters of this potential are presented in Table I. A short range repulsive term in 
the AA potential is included to simulate Pauli principle between the valence A particles. 
The core-A potential is obtained by folding phenomenological A-nucleon potential into 
the density distribution of the core nucleus. The chosen density distribution function of 
the core is given by 


p(r) 


Po 

l + exp(^) 


(31) 


with Tc = fm, a = 0.65/m, tq = 1.7fm (where r'c is termed the half density radius 

and a, the skin thickness). The the density constant po is determined by the condition 


/ 


p{r-)d?r = Ac, 


(32) 


where Ac is the mass of the core in units of nucleon mass. The value of a and tq are 
chosen following suggestions in the literature [10]. The phenomenological AN potential 
is given by 

VAN{r) = V'oexp(-rVp^) (33) 

with Vo and p adjusted to reproduce the experimental binding energy (5 a) of single-A 
hypernuclei (see Table II) in the core-A subsystem. Then the core-A potential is given 
by 

VcA{r) = J p(ri)I4v(|n - r\)d^ri. (34) 

A two dimensional plot of the AA, effective AN and core-A folded potential are shown 
in Fig. 2. The lowest eigen values (ccq) of the three body effective potential is plotted 
against the hyper-radial distance p as shown in Fig. 3. 


The strength of AN potential is expected to be weakened with the increase in mass of 
the core due to the screening or shielding effect by neighboring nucleons within the core 
when the interacting nucleon is embedded in the core. The 7r-mesic decay of A hyperon 
(A —)■ -|- tt) is predominant in the free space but tends to be suppressed in hypernucleus 
by the Pauli-exclusion principle and instead nonmesic weak process (A -|- iV N + N) 
becomes dominant with increasing mass number [7,47-51]. Thus we actually get an effec¬ 
tive AN interaction by the folding process. The parameters of this effective AN potential 
are listed in Table II. The variation of the depth of the effective AN interaction against 
the core mass (Ac) is displayed in Fig. 4. 


The truncated set of CDE for a given K^ax takes the form 




_ _ 

2/i 1 dp^ 

E ^max \ ' 

K'=0,2,.... 


n'i 


E 


5, 


Klx.^Ls{p) 


(allowed) ^ Klxi I ^l) 


KV, > Uk'uMp) = 0 


(35) 


(allowed 2, ... only for S=0, L=0). Note that the subscripts ly^ (=Gi) or I'y^ (=4i) 

have been suppressed for brevity. Eq. (12) is solved to get the ground state energy 















Figure 5: Plot of ground and excited state wavefunctions for ®^He 


F^(o) and wave functions by the hyperspherical adiabatic expansion (HAA). For 

the single-A hypernuclei, we solve Eq. (14) with V{p) replaced by Ka(p) plus centrifugal 
potential. The results are presented in the last two columns of Table II. For double- 
A hypernuclei, we use the same Ka and solve Eq. (12) for the ground state. Then 
following the prescription of sec. II, we obtain the energy and wavefunction of the higher 
excited states. The ground and excited state wave functions for aa^g, AA^^a, and aa^^ 
are respectively shown in Fig. 5, 6 and 7. The calculated values of two-A separation 
energy at Kmax = 20 for the ground state (ie. 0+), different excited states (ie. 

Om 0^) O 3 ) Gfc) and AA bond energy (AT^aa) for the ground state of aa^g, IaBg, IaHg, 
aaS, AA^^a, aa^’^ are presented in Table III. The variation of the single-A separation energy 
i?A with the mass of the core Ac is shown in Fig. 8. On the same graph, the variation of 
the two-A separation energy B\ with the mass A is also plotted. The calculated two-A 
separation energy (I^aa) and AA bond energy (AI^aa) of aa^g are in agreement with the 
experimental value 7.25±0.19 MeV [7] and 1.01 MeV [7] respectively. The predicted two- 
A separation energies for the ground and excited states of all the six doubleA hypernuclei 
considered here are depicted in the energy level diagram (see Fig. 9). 

After obtaining the three body wave function by the HHE approach, root mean square 
(rms) matter radius {Rm) of the three-body system has been calculated following the 
equation 


Rm — 


AcKc + niA < 


13 


' 12 


Ac + 2m\ 


1/2 


(36) 


whereAc, m\ are the masses of the core and the A hyperon (in units of nucleon mass) 
and Rc is the matter radius of the core determined by the relation Rc = tqA^R with 
ro = 1.70/m. The rms core-A separation is dehned as 


Rch 


< ^^13 + rU. > 


-,1/2 


12 


2 


(37) 
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Figure 10: Variation of ground state rms matter radius Rm (fm) with mass A (in nucleon 
mass unit) of the double-A hypernuclei 


The expectation value of the observable < + r ^2 > is given by the expression 


<r^3 + r^2> - T.KK'I^^LS P^dpUKi,^Ls{p)UK'i^^LS 


r/2 


X 


2a 


V cos^<h + 


IT 
( 23)1 


cos^<Fsm^<h(i<h. 


The rms separation between the valence A-hyperons (-Raa) is given by the expression 

, 1/2 


-Raa — 


<^^23 > 


(2)(<|))(2) 

(38) 

in 

(39) 


where 


<’^i3> = 


X 



P^dpUKl LS (p) Uk'I^^ ls (p) 
0 ^ ^ 

(<1>) cos'* 


(40) 


The result of these calculations are recorded in Table IV. The variation of the ground 
state rms matter radius (Rm) with mass number (A) of the double-A hyper nuclei is 
displayed in Fig. 10. 

Finally we have calculated the contribution of various orbital angular momenta at 
Kmax = 20 to the ground state probability density distribution to identify the dominating 
component of the partial waves, and the results are presented in Table V. 


IV Summary and Conclusion 

Since hyperons and nucleons both have three quark (qqq) structure (eg. p —)■ uud, 
n —)■ udd, A° —)■ uds etc.) interactions among hyperons as well as with nucleons should 


($) 
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give valuable inputs to the knowledge of strong interactions. But not much attention has 
so far been directed to the study of hyperon-hyperon and hyperon-nucleon interaction 
through the critical investigation of hypernuclei. We have undertaken a systematic study 
of the bound state properties of few light and medium mass hypernuclei to gather some 
idea about nature of hyperon-hyperon and hyperon-nucleon interactions and the internal 
structure of the three-body systems. The hyperspherical harmonics expansion (HHE) 
method adopted here is an essentially exact method, where calculations can be carried 
out up to any desired precision by gradually increasing the expansion basis. For aaHc, 
the calculated two-A separation energy i?AA and AA bond energy AI^aa (see Table III) 
both at Kmax = 20 are well within the experimental error limits 7.25 ±0.19MeF [7] and 
l.OlMeF [7]. Both values are less than onr previous calculation [52] done by nsing old 
data of [5]. The strength of the core-A potential for a hxed range parameter ( 77 ), is fonnd 
to become weaker as the mass of the core increases. And this trend is fonnd to £t a 
logarithmic fnnction of the form 

Voi^c) = loge(^c) + ^2 (41) 

with (5i = 1.55033AfeF and 62 = —43.8843MeI7. In the low mass region vq decreases 
faster than in the mid mass range before achieving satnration in the heavier mass region. 
This phenomenon may be viewed as the effect of screening of the interacting nncleon 
embedded in the core, by the surronnding nncleons. And this effect of screening increases 
as the nnmber of neighboring nncleons increases. One may expect snch a result from 
physical intuition. Following the above empirical formnla one can estimate the depth of 
the effective AN potential for the single-A hypernnclei in any mass range. The single A 
separation energy is fonnd to increase with the increase in core mass approximately 
following the eqnation 

Ba{Ac) = 6-^loge{Ac) + ^3 (42) 

with ^3 = 6.21483MeF and ^4 = —5.55812AfeF. Similarly, the trend of variation of the 
two-A separation is fonnd to follow the eqnation 

BAK{A) = 5^loge{A) + 5^ (43) 

with = 12.3892MeF and (Jg = — 11.2555MeF. Thus Fq. (42) and (43) may respec¬ 
tively be nseful in determining Ba and Baa for any one- and two A hypernnclei. Also 
the approximate size of any donble A hypernnclei of mass number (A) can be estimated 
by onr empirical formula 

Rm{A) = 5,A^^ (44) 

with ^7 = 1.25244/m and 5^ = 0.300609/m. It can be seen from energy level diagram 
(see Fig. 9) that the energy levels for the exotic isobars AA®He and AA®Be are almost 
identical. However the slight difference in their one- and two-A separation energies (see 
Colnmn 4 of Table II and III) is dne to the relatively stronger Coulomb repulsion in the 
core of the latter. The calculated values of the partial probability (see Table V) shows 
that larger contribution to the gronnd state total wave fnnction comes from the = 0 
partial wave and the rest cones from the higher even components. The contribntion 
from the odd values are zero. From the calculated values of the single- and donble- A 
separation energies (see Table II & III) and the rms matter radii (see Table IV) it may 
be conclnded that the addition of one- or two-A particle to the core makes the system 
strongly bonnd and more compact. This is qnite trne as the weakly bonnd ^He or nn- 
bonnd ^Be core becomes strongly bound on addition of one or two A hyeperon(s). Hence, 
it is qnite possible that onr predictions for the binding energies of the gronnd and excited 
states of experimentally unobserved donble A hypernnclei are valid. Thus we may say 
that that the AN force is much stronger and mnch attractive than the NN force. Perhaps 
this may be dne to the presence of the strange qnark in the A hyperon quark constituents. 
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V Figure Caption 

1. Fig. 1 Choice of Jacobi coordinates for the partition 1. 

2. Fig. 2 Plot of: (i) AA potential; (iia) effective AN potential for (iib) effective 

AN potential for ®\Zr; (iiia) core-A folded potential for ® (iiib) core-A folded 
potential for (v is in MeV and r is the relative separation) 

3. Fig. 3 Plot of lowest eigen potential uo{MeV) against hyperradial distance p for 

double-A hypernuclei (a) aa^g, (b) aa^^^ (c) AA^r 

4. Fig. 4 Plot of core-A potential depth ncA(MeV) against mass of core (Ac) 

5. Fig. 5 Plot of ground and excited state wavefunctions for aaHg 

6. Fig. 6 Plot of ground and excited state wavefunctions for aa^^^ 

7. Fig. 7 Plot of ground and excited state wavefunctions for ^AZr 

8. Fig. 8 Plot of ground state- (i) Single-A separation energy Bj^(MeV) and (ii) 
Two-A separation energy Raa (MeV), against mass of the core Ac 

9. Fig. 9 Energy Level diagram of Double-A Hypernuclei [Raa is in MeV] 

10. Fig. 10 Variation of ground state rms matter radius Rm (fm) with mass A (in 
nucleon mass unit) of the doule-A hypernuclei 
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Table II. Parameters of the AN potential and corresponding A separation 

energy in the core-A subsystems. 


System 

A-N potential 
parameters 

B^iMeV) 

Ground state 

Vo{MeV} 

Xifm) 

Expt. 

Calc. 


-43.6102 

1.41 

3.12 ±0.02 12,23 

3.1200 

;He 

-39.6328 

1.41 

5.23 ±0.00 

22 


5.2300 

;Be 

-39.3650 

1.41 

5.16 ±0.08 

23 


5.1600 

SUc 

A ^ 

-38.0436 

1.41 

17.96 ±0.00 

24 


17.9600 


-37.8497 

1.41 

19.24 ±0.00 

24 


19.2400 

-Zr 

-36.3328 

1.41 

22.10 ±0.30 

25 


22.1000 


Table III. Two-A separation energy in the ground and excited 
states of different double-A hypernuclei at Kmax = 20. 


Hypernuclei 

Bound State 

BAAiMeV) 

AHaa( 

[MeV) 


.r 

Expt. 

Calc. 

Expt. 

Calc. 

^He 

0+ 

7.25±0.19[4] 

7.2501 

0.8613 

1.01 [4] 

1.0101 

Lite 

0+ 

- 

11.5815 

- 

1.1215 

0^ 

- 

2.7056 

- 



0+ 

- 

11.4344 

- 

1.1144 


- 

2.6451 

- 


AA^ 

0+ 

- 

36.7957 

- 

0.8757 

0^ 

- 

24.4038 

- 



0^ 

- 

14.1298 

- 



03 + 

- 

07.9952 

- 



0| 

- 

02.1448 

- 


Pq 

AaL'S 

0+ 

- 

39.2357 

- 

0.7557 

0^ 

- 

27.5745 

- 



0^ 

- 

17.0814 

- 



03 ^ 

- 

10.2213 

- 



0| 

- 

04.3325 

- 


7r 

AA^r 

0+ 

- 

44.5463 

- 

0.3463 


- 

36.1553 

- 



0^ 

- 

27.0817 

- 



03 + 

- 

18.6847 

- 



Ot 

- 

11.4110 

- 



0^ 

- 

03.0006 

- 
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Table IV. The r.m.s. matter radii of the ground and excited states of 
different double-A hypernuclei at Kmax = 20. 


Hypernuclei 

Bound State 

r 

RA{fm) 

Rc\{fm) 

Rwifm) 

R{\A)c{fm) 

^AHe 

0+ 

2.3842 

3.1804 

4.1437 

2.4130 


0^ 

2.7823 

3.9536 

5.2301 

2.9652 

Life 

0+ 

2.2939 

2.9081 

3.8697 

2.1711 

0^ 

2.4952 

3.4427 

4.6361 

2.5453 

8 rL 

0+ 

2.2970 

2.9167 

3.8810 

2.1775 


0^ 

2.5005 

3.4560 

4.6542 

2.5552 

AA'J 

0+ 

3.4592 

2.9783 

4.0864 

2.1670 

0^ 

3.4666 

3.0998 

4.2331 

2.2647 


02 + 

3.4840 

3.3706 

4.6237 

2.4528 


03 + 

3.5095 

3.7331 

5.1750 

2.6909 


0| 

3.5377 

4.0998 

5.7564 

2.9196 

4^ Pq 

AAC^a 

0+ 

3.7286 

3.1140 

4.2798 

2.2623 

0^ 

3.7332 

3.2124 

4.3931 

2.3441 


02 + 

3.7450 

3.4484 

4.7293 

2.5099 


03 + 

3.7628 

3.7779 

5.2220 

2.7305 


0| 

3.7830 

4.1209 

5.7556 

2.9496 

7r 

AA^^ 

0+ 

4.9045 

3.8392 

5.3125 

2.7720 

0^ 

4.9049 

3.8581 

5.3148 

2.7969 


02 + 

4.9078 

4.0007 

5.5133 

2.8994 


03 + 

4.9124 

4.2110 

5.8141 

3.0466 


0| 

4.9179 

4.4546 

6.1672 

3.2147 


0,t 

4.9230 

4.6693 

6.4821 

3.3612 


Table V. The contribution of the orbital angular momenta to the 
probability density in the ground state of double A-hypernuclei at K^ax = 20. 


Hypernuclei 

Partial probability for = 

0 

1 

2 

3 

4 

Life 

0.996130 

0.000000 

0.003854 

0.000000 

0.000016 

Life 

0.997786 

0.000000 

0.002210 

0.000000 

0.000003 

Lb^ 

0.997757 

0.000000 

0.002240 

0.000000 

0.000004 

AA^ 

0.999900 

0.000000 

0.000099 

0.000000 

0.000001 

Pq 

AA^-a 

0.999878 

0.000000 

0.000121 

0.000000 

0.000001 

7r 

AA^r 

0.999744 

0.000000 

0.000254 

0.000000 

0.000002 
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